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ABSTRACT 



We investigate the weak gravitational lensing effect due to the large-scale structure of the 
universe on two-point correlations of local maxima (hotspots) in the 2D sky map of the cosmic 
microwave background (CMB) anisotropy. According to the Gaussian random statistics as most 
inflationary scenarios predict, the hotspots are discretely distributed with some characteristic 
angular separations on the last scattering surface owing to oscillations of the CMB angular 
power spectrum. The weak lensing then causes pairs of hotspots which are separated with 
the characteristic scale to be observed with various separations. We found that the lensing 
fairly smoothes the oscillatory features of the two-point correlation function of hotspots. This 
indicates that the hotspots correlations can be a new statistical tool for measuring shape and 
normalization of the power spectrum of matter fluctuations from the lensing signatures. 



Subject headings: cosmology: theory - cosmic microwave background - gravitational lensing - 
large-scale structure of universe 



1. Introduction 



The temperature anisotropy of the cosmic microwave background (CMB) is the most powerful probe 
of our universe because the anisotropics encode a wealth of cosmological information and the physical 
processes involved in anisotropy formation are well understood by linear perturbation theory (e.g. 



Hu 



Sugiyama & Silk 1997). Therefore, high- precision measurements of the primary CMB anisotropics enable 
us to determine a lot of fundamental cosmological parameters with unprecedented precision (Jungman et 
al. 1996| ; [Bond, Efstathiou fc Tegmark 1997| ; [Zaldarriaga, Spergcl k Seljak 1997] ). On the other hand, 
gravitational lensing due to the large-scale structure between the last scattering surface (LSS) and us 
deforms the 2D sky map of the CMB anisotropy. The weak lensing will be in the near future an effective 
tool for mapping the inhomogeneous distributions of dark matter because it is not affected by uncertainties 
in the "biasing" problem regarding the extent to which luminous objects trace the mass distribution. 
Furthermore, a great advantage in studying the lensing effect on the CMB anisotropics is that statistical 
properties of the unlensed CMB field can be fully specified once we assume the field obeys the Gaussian 
random theory ( Bardeen et al. 1986 , hereafter BBKS; Bond fc Efstathiou 1987 hereafter BE). Most 
inflationary scenarios indeed support this assumption. However, most previous studies have focused on 
theoretical predictions of lensed shapes of the two-point correlation function C{9) or equivalently the 
angular power spectrum C; of the temperature fluctuations field (Blanchard & Schneider 1987; Cole & 
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Efstathiou 1989; Sasaki 1989; Tomita & Watanabe 1989; Cayon, Martnez-Gonzalez & Sanz 1993a, b; 3eljak 
1996 ; Metcalf fc Silk 1997 ; Zaldarriaga & Seljak 1998), while some studies have shown that the lensing 



induces non-Gaussian signatures (Bernardeau 1997; Zaldarriaga 1999| ). It is therefore worth exploring an 
other statistical CMB estimator sensitive to the lensing signatures using more directly the advantage of 
predictability of the 2D CMB sky map. 

Recently, Heavens and Sheth (1999, hereafter HS) investigated the two-point correlation function 
of local maxima (hotspots) above a certain threshold in the 2D CMB map, say ^pk-pk(^)j which can be 
accurately calculated based on the Gaussian random theory once C/ is given. Their purpose was to propose 
that Cpk-pk(^') can be a powerful test of the Gaussian hypothesis of the temperature fluctuations, and 
they showed that the shape of (^pk-pkl^*) is largely different from that of C{0). ^pk-pk{^) has oscillatory 
features such as the depression around the angular scale corresponding to the 1st Doppler peak of Ci and 
a prominent peak with a damping tail toward the smaller scales than the peak-scale. The reasons can 
be explained as follows. In the derivation of i^pk-pk(^)j we need statistical properties of the gradient and 
second derivative flelds of the temperature fluctuations field in order to identify the local maxima in the 2D 
CMB map. The shape of ^pk-pk(6') thus reflects more strongly the oscillatory features of Ci such as a series 
of Doppler peaks than C{d) does, because the gradient and the second derivative fields are more sensitive 
to the oscillations of C;. The oscillations of Cpk-pkl^*) physically mean that the hotspots are discretely 
distributed with some characteristic angular separations on the LSS. Then, the weak lensing due to the 
large-scale structure causes the pairs of hotspots separated with the characteristic scale on the LSS to be 
observed with various separations in the random line of sight. This leads to the expectation that ^pk-pkl^*) 
enables us to see more directly the lensing dispersion of the specific angular separation than using C{9). 
The purpose of this Letter is therefore to develop the formula for calculating the lensing effect on ^pk-pkC^*) 



using the power spectrum approach to compute the dispersion of the lensing deflection angle (3eljak 1996) 
and to present some results in the standard cold dark matter (SCDM) model. 



2. Weak Lensing EfTect on the Hotspots Correlation Function 



During the propagation from the LSS to us, a CMB photon is randomly deflected by the inhomogeneous 
matter distributions in the intervening large-scale structure. Then, two CMB photons observed with a 
certain angular separation 9 have a different separation when emitted on the LSS. If we assume that the 
lensing fluctuations of the relative angular separation, 50{— SOi — 592), obey the Gaussian random statistics 
( ISeljak 1996| ), the ensemble average of the following characteristic function can be calculated as 



(exp[ii • {SOi - 5e2)])^Q^_Q^^^g ~ exp 



(1) 



where 56i and 862 are the angular excursions of the two photons, and ( )g denotes the average performed 
over all pairs with a fixed observed angular separation 9. The ctglI^) defined by 
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dxa-^P5{k,T)W\x,X..c){l~Mxke)), 
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where r is a conformal time, x = tq — t, Jn{x) is the Bessel function of order n, and the subscript and 
"rec" denote values at present and a recombination time, respectively. Ps{k,T) is the power spectrum of 
the matter fiuctuations, Hq = 100ft- km s^^ Mpc~^ and fimo denote the present Hubble parameter and the 
present energy density of matter, respectively. The projection operator W^XiXrcc) on the celestial sphere is 
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given by Xroc) = 1 ^ x/Xrcc in a- flat universe. The approximation used in the derivation of equation 

(0) is vahd as long as the (Jqi^{9)/9 <C 1 is satisfied, and we numerically confirmed (tgi^/O < 0.25 on the 
relevant angular scales for all normalizations of matter power spectrum considered in this Letter. 

The unlensed C{9) in the context of the small angle approximation developed by BE is given by 
Ci0) = (A(0i)A(02)) - j{ldl/2T:)S{l)MW), where |6)i - = ^ and A(0) = {T{e) - Tcmb)/Tcmb. The 
2D power spectrum of the temperature fluctuations field, S{1), is identical to Ci in the limit oil':^ 1. 

First, we briefly present the derivation of the unlensed two-point correlation function of hotspots in 
the 2D Gaussian random CMB field developed by HS. The moments of power spectrum are defined by 
erf = /g°°(WZ/27r)S'(Z)P*. At a hotspot, the gradients of temperature fluctuations field, = dA/d6^, vanish 
and the eigenvalues of the second derivative matrix, Aij = d'^A/dO^dO^, are negative. We thus need 6 
independent variables v — (A, A^, Aj,, A^x, A^y, Ayy) to specify one local hotspot. The probability density 
function of the 12 variables v = (f(i),V(2)) = {v{9i),v{92)) for the two hotspots separated with the angle 
6»(= |0i - 02 1) is defined by 



P2(t'(l),t'(2)) = 



1 



(27r)6|Af|i/2 



exp 



(3) 



by = A/ao,r]i 



{{vi — {vi)){vj — (vj))), M = det(My), and M^-^ is the inverse 
{vi) = in the present case. It is then convenient to introduce the variables defined 
Ai/(Tx,X = ~{Axx + Ayy)/a2,Y = {Axx - Ayy)/(T2, aud Z = 2Axy/(J2- AU non-zero 



where Mij is the covariance matrix: Mij 
of My . Note that 



components of Mij for two hotspots characterized by i'(i) = v{6i) and 1/(2) = v{92) and so on were explicitly 
calculated by HS. Employing the BBKS formalism, we can obtain the following unlensed correlation 
function of two hotspots above certain thresholds vi and 1^2, respectively: 
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_(iZ(2) (A(;^) - Y(j) - 2'(j-|)(A(2) - ^(2) - ^■(2)) 

) 

X P2(j^(i),^(l),y(l),2'(i),7y(i), = 0,J/(2),A:(2), Y(2),^(2),7?(2)-* =0), (4) 
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where 9^ — \/2ci/o'2, npi^{0) and npk(> i') are the number density filed and the mean number density of 
hotspots above height v, respectively. To obtain ^pk-pk(^), we performed a 6D numerical integration of 
equation (||) because the two can be done analytically. 

Next, we derive the lensed correlation function of hotspots. In this Letter, we focus on the lensing 
contributions to ^pk-pk(^) as a secondary effect which comes from a mapping of the intrinsic discrete 
distributions of hotspots on the LSS. The procedure is appropriate because a power of anisotropics generated 
by the lensing is small (Seljak & Zaldarriaga 1999; Zaldarriaga 1999| ). For this purpose, we define the 
fluctuations field of npk(0) as 5n{6) = (npk(0) — 'T-pk)/^pk- Likewise, since the gravitational lensing makes 
us to observe, at a certain angular position 0, a hotspot at a different position 6 + 56 on the LSS, the lensed 
(observed) fluctuations field dh{9) can be expressed as Sh{9) = 6n{6 + 66) = j d^l / 5n{l) exp[ji- (0-1-50)], 
where 5n{l) is the Fourier component of dn{9). If assuming that the lensing fluctuations and the intrinsic 
CMB anisotropics field are statistically independent, then the lensed two-point correlation function of 
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hotspots can be calculated with help of equation (||) as 

1 



(<5n(0i)<5n(02)) 



where Io{x) is the modified zeroth-order Bessel function 



y"d^'0'epk-pk(e')exp 



+ 1 



(5) 



This equation indicates that the lensing 
contribution at a certain scale arises only from the lensing dispersion (Tgl(^) at the same scale alone. 
Hence detecting the scale dependence of the lensing contributions to ^pk-pk(^) enables us to reconstruct the 
lensing dispersion, more importantly the projected matter power spectrum by equation 



3. Results 

In the following, we show some results in the SCDM model with fimo = 1- The Hubble parameter and 
the present baryon density are h — 0.6 and Obo^-^ = 0.015, respectively. To compute the intrinsic CMB 
angular power spectrum, S{1), we used the CMBFAST code developed by Seljak & Zaldarriaga (1996). As 
for the matter power spectrum in the lensing dispersion (|2|), we employed the Harrison-Zel'dovich spectrum 
and the BBKS transfer function with the shape parameter in Sugiyama (1995). A free parameter in our 
model is only the normalization of the present-day matter power spectrum which is conventionally expressed 
in terms of the rms mass fluctuations within a sphere of 8ft-~^Mpc, i.e., as- 

In Fig. ^ we show the two-pint correlation function of hotspots of height above the threshold 
V = \ with and without gravitational lensing effect, ^^^pk(^) and ^pk-pk(^'), respectively. Note that the 
intrinsic ^pk-pk(^) can be calculated by using only Ci. The mean number density of hotspots is then 
"-pk(> v) ~ 1.16 X lO^'rad"^. The dependence of the lensing effect on the normalization of mass fluctuations 
is demonstrated by the four choices of o-g — 0.5,1,1.5,2, where a% = 0.5 and 1 roughly correspond to 
normalizations of the cluster abundance (Eke, Cole & Frenk 1996; Kitayama fc Suto 1997] ) and the COBE 



measurements (Bunn & White 1997), respectively. The figure clearly shows that ^pk-pk(^) has more 
oscillatory features than C{9). This is because the gradient and second derivative fields of A(0) used in 
the derivation of ^pk-pk(^) have power spectra of CiI'^/{2tt) and Cil^/{2n) per logarithmic interval in I, 
respectively, and thus they strongly enhance the oscillations of PCi/{2tt), while C{9) is fully determined 
by PCi/{2tt). If using the relation / « tt/O, there are indeed correspondences between both oscillations of 
Cpk-pk(^) and l^Ci such as the depression at 6* w 75' corresponding to the scales around the 1st Doppler 
peak of Ci, a prominent peak at ~ 13', and a damping tail ai 9 < 10' associated with the Silk damping. 
Moreover, Fig. |^ shows that the gravitational lensing appears as a larger smoothing at scales where 
Cpk-pk(^) fias more oscillatory features, since pairs of hotspots separated with the characteristic angular 
scales corresponding to the oscillations of ^^pk-pk are redistributed by the weak lensing and thus power of 
Cpk-pk(^) is transfered from (or to) the scale to (or from) nearby scales. Especially, it should be stressed 
that even the lensing effect on the depression at 9 ^ 75' associated with the 1st Doppler peak of C; is 
distinguishable where the deviation S£_ (= (■^^k^pk ~ fpk-pk)/'?pk-pk) reaches up to « 20% for erg = 2.0 
while SC < 1.0% in all cases (see Fig. ^). Since both C^^pk and Cpk-pk give the same results as the 
shape becomes fiat on larger scales {9 > 80'), the large scale amplitude of Cpk-pk gives the normalization 
independent of the lensing contributions. We also demonstrate the nonlinear corrections to the matter 
power spectrum to clarify how the nonlinear evolution of matter fluctuations affects our analysis. Using 
the semi-analytic formulae developed by Peacock & Dodds (1996) for the extreme case of as = 2, where 
the universe has too large power of the present-day small-scale matter fluctuations, we found that the 
nonlinear correction is small at 6* > 5'. Most importantly. Fig. |l] indicates that the magnitude of the lensing 
contributions to ^pk-pk is fairly sensitive to the amplitude of erg- Furthermore, measuring scale-dependence 
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of the Icnsing contributions such as S^{d) a,t 9 = 2', 16', 36', 44', 75' will allow us to reconstruct the shape of 
projected power spectrum of matter fluctuations in the range of 3/i~^Mpc ^ A ^ 50/i~^Mpc. 



4. Discussions and Conclusions 

In this Letter, we have investigated the weak lensing effect on the two-point correlation of hotspots 
in the CMB. One of the great advantages in studying the lensing effect on the CMB anisotropics is that 
statistical properties of unlensed CMB field are fully predictable by the Gaussian random theory. For 
example, Bernardeau (1998) have investigated how the lensing changes probability distributions of the 
ellipticity defined from the local temperature curvature matrix, and shown that the method can be a 
statistical indicator sensitive to the lensing signatures. 

It is now expected that the weak lensing can be a powerful probe for measuring the power spectrum of 
matter fiuctuations. In case of the normalization based on measurements of the large angular scale CMB 



anisotropies such as COBE normalization (e.g., Bunn & White 1997), we need to employ the very accurate 



matter transfer function which highly depends on a lot of cosmological parameters in order to determine 
the amplitude of small scale matter fiuctuations. In fact, even most sensitive future experiments such as 



the Planck Surveyor (Bersanelli et al. 1996) will constrain as up to only 20% accuracy (see, e.g.. Bond, 
Efstathiou fc Tegmark 1997). We thus wish to propose that ^pk-pk(^') can be a new effective statistical tool 



for measuring amplitude and shape of the matter power spectrum. 

For this purpose, we have to perform quantitative predictions whether or not the future satellite 



missions, MAP ( Bennet et al. 1996 ) and Planck, can measure the lensing contributions to ^pk-pk(^) and 
constrain erg with high accuracy. In particular, we have to estimate not only uncertainties by the cosmic 
variance but also observational errors in identifying angular positions of hotspots in a realistic 2D sky map 
of the CMB anisotropies. It is then important to take into account the following experimental effects. 
(1) The finite beam size of detectors may incorporate some hotspots within the beam size resulting from 
smoothing of the lensed CMB anisotropies filed. The effect hence decreases the mean number density of 
hotspots and suppresses two-point correlations of hotspots on scales below the beam size. (2) The detector 
noise may make spurious hotspots in the observed CMB map. HS estimated the errors for measurements 
of the unlensed hotspots correlation functions expected from MAP and Planck by using numerical sky 
realizations, and indeed showed that the theoretical predictions including the above experimental effects 
are in good agreement with the numerical results. Following their method, we are now investigating the 
detectability of the lensing signatures to ^pk-pk(^), and have obtained the preliminary result as follows. 
Supposing the data produced by the Planck specification with 65% sky coverage, the expected signal to 
noise ratios for the lensing signatures are S/N w 1.3, 1.8, 1.5 and 3.1 at 6* = 16', 36', 44' and 75', respectively, 
for (78 — 1.5 model. Furthermore, we expect that the statistical significance of S/N could be increased by 
combining measurements of the two-point correlations of coldspots. The observed correlations of coldspots 
should be identical to that of hotspots since they are subject to the Gaussian random statistics. These 
works are now in progress and will be presented in detail elsewhere. 

We finally comment on possibilities that the lensing effect on ^pk-pk breaks the so-called geometrical 



degeneracy ( Bond, Efstathiou fc Tegmark 1997 ; Metcalf fc Silk 1998 ; gtompor fc Efstathiou 1999 ) inherent 



in cosmological parameter determinations from measurements of the CMB anisotropies. We have confirmed 
that the global shape of ^pk-pk(^) is sensitive to fundamental cosmological parameters such as h, fib, ^mo 
and so on. We thus expect that detecting 6'-dependence of the lensing contributions to C^k^pk(^) '^^^ yield 
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an additional constraint on the erg — Qrno plane if we assume the shape of matter power spectrum. Hence, 
the degeneracy could be broken by observations of the CMB anisotropics alone, without help of other 
astronomical measurements. 
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Fig. 1. — The correlation function of hotspots vs. 9 with and without the gravitational lensing effect 
in SCDM model. The solid line is the unlensed correlation function, while dotted, short-dashed, long- 
dashed and dotted-dashed lines represent the corresponding lensing cases with as = 0.5, 1.0, 1.5 and 2.0, 
respectively. The long-dotted-dashed line denotes the result including nonlinear correction with erg = 2.0 
(see text), h = 0.6 and f^bo^^ = 0.015 are assumed. The threshold of hotspots is = 1. The mean density 
of hotspots is then fipk = 1.16 x lO^'rad"^. 
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Fig. 2. — The two-point correlation functions of hotspots and the temperature anisotropics field around the 
fst Doppler peak with and without the gravitational lensing effect. 



